Quantum gates using electronic and nuclear spins of Yb+ in a magnetic field gradient 
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An efficient scheme is proposed to carry out gate operations on an array of trapped Yb"*" ions, 
based on a previous proposal using both electronic and nuclear degrees of freedom in a magnetic 
field gradient. For this purpose we consider the Paschen-Back regime (strong magnetic field) and 
employ a high-field approximation in this treatment. We show the possibility to suppress the 
unwanted coupling between the electron spins by appropriately swapping states between electronic 
and nuclear spins. The feasibility of generating the required high magnetic field is discussed. 
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I. INTRODUCTION 

Quantum computing (QC) based on nuclear spins has 
attracted considerable interest over the past decades due 
to the long coherence time of nuclear spins. Since they 
only weakly interact with their environment, nuclear 
spins are well suited for storing quantum information, 
and, for the same reason, difficult to manipulate. As a 
result, to carry out QC, one may work with ensembles 
of nuclear spins [l| or employ a hyperfine interaction to 
manipulate individual nuclear spins by electron spins . 
The former, using mature techniques of nuclear magnetic 
resonance, has become a test bed for models and schemes 
of QC and quantum simulations 0, 0] , while the exploita- 
tion of the latter is still at an early stage. 

Furthermore, there have been proposals combining nu- 
clear and electronic spins in solid-state systems, such as 
doped silicon substrates [5| and doped fullerenes @. To 
distinguish nuclear and electronic degrees of freedom, one 
has to introduce a strong magnetic field and work in the 
Paschen-Back regime. An experimental demonstration of 
the manipulation of a nuclear spin ensemble with about 
10,000 ions in a Penning trap under magnetic field 0.8 T 
was reported by Bollinger et al. Q. The manipulation 
of individual nuclear spins had only been achieved in di- 
amond nitrogen- vacancy centers, with qubits encoded in 
the ^'^C or ^^N nuclear spins near the electron spin Q. 

Here, we propose a scheme, using trapped ^^^Yb+ ions, 
to encode qubits in both electronic and nuclear spins 
of trapped atomic ions for QC and quantum simula- 
tions, following a previous idea Different from al- 
ready accomplished experimental work with trapped ions 
[loj . we consider quantum logic operations on ions in 
the Paschen-Back regime with qubits encoded in nuclear 
spins / and auxiliary qubits in electron spins S = 1/2. 
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This combines the long decoherence time of nuclear spins 
with efficient manipulation and readout using electron 
spins. Quantum information is stored in nuclear spins 
and is only swapped into electronic spins for single-qubit 
gates and conditional quantum dynamics with two and 
more ions. Thus, quantum information remains well pro- 
tected from ambient noise fields that otherwise would 
give rise to decoherence. 

Swapping quantum information between nuclear and 
electron spins is accomplished using microwave radiation. 
Subsequent conditional quantum dynamics between elec- 
tron spins and individual addressing of electron spins may 
be done using laser light [lllll^l, or, in the presence of a 
spatially var ying magnetic field, using again microwave 
radiation [13l - [l9i. A magnetic field gradient induces spin- 
spin coupling between electronic spins that can 
be used for quantum logic gates, and, in addition allows 
for ions to be addressed in frequency space. The latter 
approach is useful, since it avoids technical and funda- 
mental difficulties when using laser light for coherent op- 
erations [TgI . [isl . [20j . In addition, it allows for condi- 
tional quantum dynamics without stringent re quir ements 
on the cooling of the ions' vibrational motion |21| . 

The key point of our proposal is the suppression of 
unwanted coupling between electronic spins in the mag- 
netic field gradient by swapping quantum information be- 
tween nuclear and electron spins by microwave radiation. 
As a result, the overhead operations, such as refocusing 
pulses, in [9| is no longer necessary. In what follows, we 
first investigate the use of nuclear and electron spins of 
atomic ions with and without the high-field approxima- 
tion, from which we know how well the quantum gate 
is performed in the real experimental situation. Also, 
a spatially varying magnetic field is included. Then we 
show how two-qubit gating is achieved without the need 
for compensating for unwanted couplings between elec- 
tron spins. Furthermore, a detailed discussion of how the 
required magnetic field and the magnetic gradient can be 
achieved is given. 

For concreteness, the present work considers as an ex- 
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FIG. 1: Energy-level schemes of ^'''^Yh+ (Left) and "iYb+ 
(Right), where the energy gaps are not drawn to scale. 



ample the use of ^^^Yb^ ions for quantum information 
science pTl |22| - [27| . ^''^Yb+ features a nuclear spin of one 
half and thus provides the simplest hyperfine structure 
with several potential qubits where the experimenter can 
choose both magnetically sensitive and insensitive qubit 
transitions (to first order). The hyperfine qubits can be 
directly manipulated by a resonant microwave field or by 
using an optical Raman transitions. Recent interest in 
trapping Yb"*" is, to some extent, motivated also by the 
fact that the experimenter benefits from low priced diode 
lasers for photoionization and for all transitions rel- 
evant to Doppler cooling, qubit initialization and state 
selective detection(See Fig. 1). 

The paper is structured in four sections. In section II 
we present the Hamiltonian describing a linear Coulomb 
crystal of ions and justify the high-field approximation. 
In section III we will demonstrate efficient QC operations 
using both nuclear and electron spins, the experimental 
feasibility of which is discussed in detail in section IV. 
We give a brief summary in the last section. 



II. THE SYSTEM AND THE HAMILTONIAN 

We consider an array of trapped ions in a linear trap in 
the presence of a magnetic field gradient, whose Hamil- 
tonian in units of ?i = 1 is written as, 

H^Y. ^sSl + + aY^siii + sill + Sill) 

i i i 

where Sk and Ik (k = x,y,z) are, respectively, the spin 
operators of the electron spin {S = 1/2) and the nu- 
clear spin /. For ^^^Yb+ we have / = 1/2, the hyperfine 
coupling constant is A = 12.645 GHz and the Larmor 
frequencies are given by ils ~ gsfJ-sB = 28B GHz and 
fij = gnfJ,BB ~ —7.5-8 MHz with B the strength of the 
magnetic field in Tcsla experienced by the ion. Jij is 



the coupling between the electron spins of the trapped 
ions under the magnetic field gradient p^| - [l7j . We leave 
out the vibrational modes here, because additional radi- 
ation fields applied to swap information between nuclear 
and electron spins only drives carrier transitions, that 
is, these fields do not couple vibrational and spin states. 
The nuclear spin couplings arc also neglected as they are 
very small compared to other terms. 

The magnetic field is applied along the trap axis, so 
the zth ion experiences the magnetic field with 

[Bo + bz%, 

with Bq the strength of the magnetic field at the origin, 
b — dB/dz the magnetic field gradient, and the unit 
vector along the trap axis. For the spin-spin coupling we 
have 
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where m is the mass of trapped ions, = v^^i with 
the axial frequency of the trap and /i; the eigenvalue of 
potential Hessian matrix. D is the unitary transforma- 
tion matrix that diagnonalizes the Hessian matrix and 
rig depends on magnetic gradient b. With respect to the 
original expression of Jij in [13, Eq. (2) seems for- 
mally larger by 4 times, which is because we use angular 
momentum operators here instead of the Pauli operators. 

For our purpose, we first consider the single-ion case 
to justify the high-field treatment. In the Paschen-Back 
regime, the Hamiltonian of a single ion is obtained by 
reducing Eq. (1), 

Ho = nsS, + Qil, + AiS.,h + Syly + SJ,). (3) 

Assuming a magnetic field _B = 1 T, we plot the energy- 
level structure determined by Eq. (3) in Fig. 2 for 
^''^Yb"'", where the eigenstates include some superposi- 
tions due to the x- and y-terms of the hyperfine interac- 
tion. In contrast, the conventional treatment, to simplify 
the problem, is the exclusion of the x- and y-tcrms of the 
hyperfine interaction under the high-field approximation, 
i.e., 



Hi = nsS, + nih + AS J, 



(4) 



Since each term in this Hamiltonian is diagonal, the 
eigenstates of the Hamiltonian are ones of Sz or 1^ . 

In order to carry out single-qubit gates and con- 
ditional quantum gates with nuclear spins (that 
arc used as a quantum memory), it is necessary 
to transfer the nuclear spin's state to the elec- 
tron spin and vice versa. Since the SWAP gate 
could be performed by appropriate CNOT gate se- 
quences, e.g., SWAP = CNOT/sCNOTs/CNOT/s = 
CNOTs/CNOT/sCNOTs/ [| with CNOTah implying 
the control a and target 6, we consider below the nec- 
essary CNOT gates, which could be accomplished by ra- 
diating the ion with appropriate n pulses The key 
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FIG. 2: (Color online) Angular momentum states in the 
ground state of a ^'^Yb^ ion in a magnetic field B =1 T 
(i.e. in the Paschen-Back regime). The levels drawn with the 
black solid lines and the red dashed lines are, respectively, 
from the exact Hamiltonian Eq. (3) and the approximated 
Hamiltonian Eq. (4). 



point is the consideration of the level shifts due to hy- 
perfine interaction. 

Taking a CNOT5/ gate as an example, under the high- 
field approximation with the magnetic field 1 T, we may 
radiate the ion by a 6.31 GHz microwave pulse, yielding 
the flip between |l/2,l/2) and |l/2,-l/2) (See energy 
levels in red by dashed lines in Fig. This pulse does 
not lead to transitions between the levels | — 1/2, 1 /2) and 
I — 1/2,-1/2) due to detuning. As a result, the flip of 
the nuclear spin is controlled by the electronic spin. The 
idea is easily extended to perform a CNOT/g with the 
nuclear spin as the control qubit. If the exact treatment 
is used instead of the high-field approximation, however, 
the involvement of x- and y-terms of hyperfine interaction 
in Eq. ^ makes the CNOT gates considered above less 
perfect. 

As shown by the black solid lines in Fig. [2l to achieve 
the transition |l/2,l/2) o |l/2, — 1/2), we may employ 
the pulse with frequency 4.95 GHz, which actually leads 
to 

|l/2, l/2)e. o 0.9776|l/2, -l/2)e..+0.2103|-l/2, 1/2)^., 



means the state under exact evolution. So 

1/2,1/2) 



where | . . 

the ion would leak to the unwanted state 
state with probability 0.04. 

To have a preferable quantum gate in a realistic ex- 
periment under the high-field approximation, we have to 
first justify the condition for Eq. (4). To this end, we 
require the level splittings under the high-field approxi- 
mation to be identical to the exact situation, which makes 
the theoretical treatment closer to the realistic operation. 
So we introduce effective gyromagnetic ratios to replace 
the natural gyromagnetic ratios in the treatment of Eq. 
(4). By numerics, we have found the effective Larmor 
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< Bq < 5 T. Effective gyromagnetic ratios 
read 7^, « (28.1 + S.Se^i-^^"/^!) and 7} « -(0.085 + 
5 5g-i.5Bo/Bi) (-^i^j^ ^ I rj,y ^^^^ ^j^ggg effective Lar- 
mor frequencies, Eq. (4) becomes 



Hi — fl'gSz 



(5) 



Eq.® is a good approximation to Eq. (3) which gives 
nearly the same energy levels. In Fig. 3 we show a high- 
ter fidelity is achieved for a CNOT gate with increasing 
magnetic field. 
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FIG. 3: (Color online) Fidelity C of the CNOTg/ gating with 
respect to the magnetic field strength B divided by hyperfine 
coupling constant A. Here C = exact)^^ with \^) and 

exact) the evolved wavefunctions under Eq. ([5]) with effec- 
tive Larmor frequencies and under Eq. ((3)1. This result with 
scaled magnetic field can be applied to other ions with 1=1/2 
and S=l/2 in a good approximation (i.e., neglecting the dif- 
ference between nuclear Larmor frequencies of different ions). 
For ^"Yb +, the magnetic field changes within the range (0.9 
T, 5 T). 

From now on we use the approximate Hamiltonian 
Eq. (|4]) with effective Larmor frequencies to simplify the 
treatment in multi-ion situation. We have to emphasize 
that the purpose of the approximation we employ is, on 
the one hand, to keep consistent with the conventional 
treatments in previous works 0, Q, on the other hand, 
to have a clear physical picture for our gate operations. 
Since the nuclear and electronic spins are never com- 
pletely decoupled in real case, we justified the approx- 
imation in the above treatment to try to find a trade-off 
for achieving accurate and coherent gate operations. 



III. QUANTUM GATING USING S-I SWAP 

Consider a string of trapped ions in the presence of a 
spatially varying magnetic field along the z-direction, as 
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described by Eq. (1). In the Paschen-Back regime, we 
may neglect the x- and y-terms in the hyperfine interac- 
tion. So the Hamiltonian is reduced to 



qi qj 



. . (6) 

The nuclear spins, due to negligible coupling with each 
other, remain the same as in the single ion case. But we 
have to pay more attention to the electron spins, which 
are coupled due to the magnetic field gradient. Because 
of these J couplings, the transition frequency of a given 
ion depends on the electron spin states of others. So 
with an increasing number of ions the spectrum of the 
ion chain becomes more complex. 

A good candidate system for gate operations should 
have the coupling between qubits well controlled. In the 
absence of a magnetic field gradient, trapped ionic qubits 
interact by coupling to the common vibrational modes 
mediated by suitably tuned laser light [Tol - [l^ . Here, in 
the presence of a magnetic field gradient, the electron 
spins' coupling, that reaches well beyond nearest neigh- 
bors, is to be used for conditional quantum dynamics. 
Other QC proposals that make use of nuclear and elec- 
tron spins usually assume only nearest-neighbor coupling 
0. For our present trapped ion model, however, the in- 
teractions between the ions are significantly beyond the 
nearest-neighbor couplings. 

Previously proposed solutions to this problem include: 
(1) refocusing operations applied simultaneously with the 
gating pulses on the trapped ions @ and (2) additional 
potentials applied on the trapped ion [l^ [23| . The for- 
mer solution is based on exact knowledge of the undesired 
coupling, and overhead in this method increases quickly 
with the number of qubits. The latter requires micro- 
structured electrodes traps to shape the effective poten- 
tial confining the ions. In order to produce sizeable J 
couplings in the system, the electrodes' axial extension 
should be of the order of 10 micrometers or smaller and 
the distance between the electrodes' surface and the ions 
needs to be of similar magnitude. 

In what follows we present an alternative approach to 
accomplish high-fidelity two-qubit gating making use of 
the two spins available in each trapped ion. Since both 
the electron spin and the nuclear spin are initially polar- 
ized, the quantum information during the QC implemen- 
tation is only stored in one of them and the other remains 
well polarized. When an ion is active, i.e., operated for 
gating or readout, the quantum information is encoded 
in the corresponding electron spin and the corresponding 
nuclear spin remains well polarized. When the ion turns 
passive, the quantum information is swapped to the nu- 
clear spin for storage, and the electronic spin becomes 
well polarized. 

Consider two active ions in an array of ions. Since the 
electronic spins of all other ions are well polarized and 
their nuclear spins, encoding the quantum states, have 
no interaction with other ions, the two active electronic 
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FIG. 4: (Color online) Spectrum of trapped ions with mag- 
netic field Bo = l T, field gradient b=500 T/m, and axial trap 
frequency 1^2=600 kHz, where (a) is for three ions and (b) 
for four ions. The lowest trace extends over all resonances; 
the middle trace shows the splitting due to a magnetic field 
gradient, and the upper trace depicts the splitting due to J 
coupling between the ions. Resonant frequencies are marked 
in lowest and middle traces, and splittings are labeled in the 
middle and upper traces (Some splittings with no labeling due 
to symmetry). The ions (Nos. 2 and 3) in red stand for active 
ions and others in gray for passive ions. The frequency lines 
in red are the ones left after the passive electronic spins are 
polarized. 



spins only experience a frequency shift by the polarized 
electronic spins of the rest passive ions. Compared to the 
case with the passive electron spins in arbitrary super- 
position, this scheme makes the spectra much simpler. 
As a result, quantum gating would be much easier since 
refocusing unwanted interactions or operations with lo- 
cally shaped electrostatic potentials in micro-traps arc no 
longer necessary. In addition, quantum states stored in 
nuclear spin degrees of freedom are more robust to de- 
coherence than in electronic counterpart, which helps to 
store quantum information in higher fidelity. 

Fig|3] demonstrates the cases with three- and four-ions 
as examples, where two of them are active and the rest 
are passive. By polarizing the passive electron spins, 
there are only two frequency lines left for the active spins, 
which arc shifted with respect to the original positions by 
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corresponding J couplings. 

We have simulated the two-qubit gating for three and 
four ^^^Yb^ trapped in a line under a magnetic field gra- 
dient, as shown in Table I. The gating time is dependent 
on the magnetic field gradient, but not on the magnetic 
field itself. A stronger magnetic field is required just for 
a better gate manipulation. 

TABLE I: Two-qubit CNOTsjSa gating time T for three and 
four trapped Yb"*" in a strong magnetic field Bq = 1 T with 
axial trap frequencies Vz = 600 kHz and 200 kHz for different 
magnetic field gradients h and distances Azmin. J represents 
the nearest-neighbor coupling for Jij in Eq. (6) and data for 4 
ions are regarding the middle two ions. N stands for number 
of ions. 



i^.(kHz) 


N 


b (T/m) 


AZmin (^m) 


J (kHz) 


T (ms) 




3 


50 


4.15 


0.0444 


70.8 
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too 


4.15 


0.178 


17.7 


600 


3 


300 


4.15 


1.60 


1.97 
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50 


3.50 


0.0368 


85.2 
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too 


3.50 


0.147 


21.3 
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300 


3.50 


1.33 


2.37 
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50 


8.63 


0.399 


7.87 
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too 


8.63 


1.60 


1.97 


200 
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300 


8.63 


14.38 


0.218 
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50 


7.28 


0.332 


9.47 




4 


100 


7.28 


1.33 


2.37 




4 


300 


7.28 


11.94 


0.263 



IV. EXPERIMENTAL FEASIBILITY 

In order to implement QC and quantum simulations 
with nuclear spins of trapped ions as described above, 
a strong and highly stable magnetic field is required. If 
conditional quantum dynamics is carried out based on 
magnetic gradient induced coupling [l^, [13, EB-fSl using 
microwave or radio frequency radiation, instead of laser 
light, then the applied magnetic field in addition needs 
to vary spatially. In addition, a magnetic field gradient 
allows for individual ion addres sing in frequency space 
using rf or microwave radiation jl3l [Tgj. Below, we will 
discuss the feasibility of creating the required strong field 
using Yb"*" ions as a concrete example. 

Qubits may be encoded in the simple hyperfinc struc- 
ture of the isotope ^^^Yb+ . Gradients can be achieved 
by using permanent magnets, for example, in quadrupole 
configuration or by using a pair of anti-Helmholtz coils or 
shaped planar current geometries known from magnetic 
traps for neutral atoms jsol . [3T| . The field noise can be 
reduced by superconductive materials [l^. 

In micro-structured traps (two- and three- 
dimensional), the required magnetic gradient extends 



only over a limited volume, and thus does not neces- 
sarily require strong fields. The optimization of the 
field geometry using current carrying micro-structures 
to reach a maximum gradient with a limited current (or 
dissipated power) appears not too different from the task 
of optimizing for maximum field. To reach magnetic 
fields in the Tesla range, however, massive cooling of the 
current carrying structures would be necessary. 

A three-dimensional ion trap has been designed that 
allows for generating gradients of up to 100 T/m j29j . 
Details on this trap will be given elsewhere. For neu- 
tral atom trapping, two-dimensional current structures 
were exploited to create flexible magnetic field configu- 
rations and magnetic gradients |30l - [33 |. A good thermal 
contact between current carrying structures and the sub- 
strate allows to efficiently remove any thermal intake due 
to ohmic heating, resulting in enormous possible current 
densities of jmax ~ 10^^ A/m^ [s^] and allows for versatile 
and fast switchable fields and gradients. With appropri- 
ate cooling of the substrate, we expect gradients in the 
range of 100-300 T/m to be possible for two-dimensional 
traps that arc currently under development in our labo- 
ratory at Siegen. 

An alternative straightforward solution for producing 
both strong magnetic fields and high magnetic gradients 
would be a pointed yoke which, in its proximity, would 
create a combination of both. It is desirable, however, to 
create homogeneous and inhomogeneous parts indepen- 
dently: the gradient is necessary for the addressing of 
single ions and the coupling between ionic qubits, but 
it impedes the efficient cooling of the ion chain as a 
whole, since the microwave transition, which is also re- 
quired during cooling to avoid optical pumping would be 
different for each ion. It can thus be advantageous to 
switch the gradient on during manipulation only, and set 
it to zero during cooling, and potentially during read-out. 
The strong offset field, however, that would indiscrimi- 
nately shift all resonance frequencies over widely spread 
frequency bands upon switching, remains on at all times. 
Therefore, in what follows, we will focus on an indepen- 
dent creation of homogeneous and inhomogeneous mag- 
netic fields. 



A. Possible approaches to creation of high 
magnetic fields 

Methods to generate strong magnetic fields suitable for 
ion trap QC include: 

• Permanent magnets produce stable and low noise 
fields, which are simple and inexpensive, but cannot 
be switched on and off nor be tuned directly - the 
field on a given point however, could be changed by 
changing the position or orientation of a magnet. 

• Current carrying structures on the other hand can 
produce time dependent fields, but require high 
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power current supplies, and in most cases stabi- 
lization and cooling. 

• Superconducting current carriers have a better in- 
herent stability but require a high initial experi- 
mental effort and costs for setting up a cryostat, 
and thereby in the long term exhibit high operat- 
ing costs. 

When using type-II super-conductors (s^j , the per- 
sistence of the magnetic field in superconducting 
coils is usually viewed as an advantage and reduces 
ac field noise, which, in terms of coherence time 
is, of course, desirable. On the other hand, per- 
manent magnets offer intrinsic low field noise and 
the persistence of superconducting magnets can be- 
come an obstacle, in case one intends to modulate 
the magnetic field periodically using additional cur- 
rent carrying structures. This can be useful, for 
instance, in order to insert temporal phases with a 
homogeneous magnetic field and thus homogeneous 
or global cooling. Slow variations are possible, but 
not on the time-scale of the typical repetition rate 
(order of 100 Hz) of data taking. 

For a static and homogeneous offset field, we will fo- 
cus on permanent magnets, where much progress has 
been made during the last decades [13]: Not only did 
the maximum remanence of commercially available per- 
manent magnets increase substantially (the remanence of 
Nd2Fei4B can reach values of up to Br = 1.22 T [H), 
but also progress was made in the task of maximizing the 
field with a given magnet material by choosing a suitable 
mounting geometry. Our investigations concerned with 
the creation of strong magnetic fields, therefore, focus on 
the usage of optimized magnet arrangements to exceed 
the surface flux of a single magnet. This can be achieved 
for example by pol e and yoke design or in Halbach ar- 
rangements [3a, |33| . 

B. The Halbach structure 

Ideally, a Halbach structure consists of an infinitely 
long magnetized cylinder of continuously varying magne- 
tization direction with inner diameter and outer diam- 
eter Tq. The magnetization of a Halbach dipole points 
along the angle 2(j) for an infinitesimal cylinder segment 
at angle cf). This results in a cancelation of fields outside 
of the cylinder and a homogeneous magnetic field inside 
the cylinder with the magnitude 

B = Br\n(^^y (7) 

with Br being the remanence. 

A conclusion for the generation of high magnetic fields 
can be seen from this idealized analytic expression Eq. ([T] 
for limited outer dimensions, a smaller inner diameter 
allows for larger field strength. The trap structure and 




FIG. 5: Schematics for the ideal Halbach configuration (left) 
and a segmented approximation with A'^ = 16 segments 
(right), which theoretically delivers more than 97% of mag- 
netic field strength in the inner cylinder as compared to the 
ideal structure. Black arrows indicate the local direction of 
the magnetic field. 

vacuum housing inside the Halbach cylinder cannot be 
made arbitrarily small, effectively limiting the achievable 
field strength. 

This structure can be replaced, for the ease of fabrica- 
tion, by a segmented cylinder, made of homogeneously 
magnetized segments, as shown in Fig. [51 where the field 
inside becomes 

and with = 16 segments one can reach already 97% of 
the field strength of the idealized case. 

The finite length of any real Halbach structure, too, 
contributes to a reduction of the B field according to 

B = Br \n(^^^ - Brfizo), (9) 

where the reduction factor /(zq) depends on the length 
zq and the radii and of the cylinder as 

We carried out numerical simulations for a structure 
with 16 segments with a cylindrical inner volume of a 
diameter 5 cm using NdFe35, with a retentivity of Br ~ 
1.23 T. A very homogeneous field of 2.3 T is obtained 
when the outer diameter is limited to 50 cm (see Fig. [6]). 

Even stronger fields can be achieved in three- 
dimensional structures which follow the same concept, 
namely Halbach spheres, at the expense of constrained 
optical access to the high-field region [l^ . The theoreti- 
cal field for such an arrangement is given by 

B^^Brln(^^y (10) 

which is already larger by a factor 1/3 than the field 
created by a comparable cylinder, but the low optical 
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FIG. 6: (Color online) Numerical simulation with a Halbach 
dipole ring with an inner diameter of 5 cm and an outer di- 
ameter of 50 cm, yielding a homogeneous offset field of 2.3 T, 
in contrast to the analytical model, yielding above 2.8 T. In 
addition, a set of flux lines is shown. 



access makes this choice less attractive, unless the whole 
detection is placed inside the sphere. 

In such structures, the field can exceed the maximum 
coercivity of the permanent magnet and locally reverse 
its magnetization (if aligned unfavorably) , thus imposing 
another practical limit on attainable field strengths. This 
can be avoided, by replacing parts of the magnets with 
materials with high coercivity (and often lower rema- 
nence), and in this way, magnetic fields exceeding 4.5 T 
have been created [sll . 

The maximum attainable field strength is limited by 
material properties as remanence and coercivity, which 
are usually functions of temperature. For example for 
NdFeB, the temperature coefficient of the remanence is 
-0.1 % K~^, and the temperature coefficient for the co- 
ercivity is -0.6 % K~^, allowing for substantial improve- 
ments even if cooling only to liquid nitrogen temperatures 

MM. 



The Hamiltonian reads H = T-Lq + Hi coafltt, with 



Ho 



_ (p-eAf 



and 



Hi = -\mn^M^^ - y^) , 



(11) 



where a and q are the usual stability parameters charac- 
terizing the trapping potential [4^ . 

The corresponding Schrodinger equation can be solved 
with the ansatz 



*(x,2/,t) = $(.T,y,i)e-*"(*)^'("^^^^ 



(12) 



of a slowly varying wave function $(x, y,t) and the 
quickly oscillating phase a = sin ^Itt- Averaging the 
time-dependent Schrodinger equation over the interval 
St = 2TTn^^ and taking <i> to be constant over this period 
yields 



ih— = H$ 

at 



(13) 



with 



n = —ipi+pt,) + -mutix^+y') + -u,{p,y~pyx) (14) 
where the following averages have been used: 



adt 



ik [^""^t ' ^2 

27r Jo 



a^dt 



and 



2tt Jo 
VHi = - 

Introducing the standard creation and annihilation op- 
erators 



and 



yields 



x,y 



% — hwr{a\cLx + a\ay + 1) + -huJc{a\ay — a^Gx) , (15) 



C. Effective Potential with magnetic field 

Given the substantial strength of magnetic fields con- 
sidered here, its impact on the ion's motion is a priory not 
necessarily negligible. To assess the impact of a strong 
magnetic field, we consider here the dynamics of an ion 
with mass m and charge e in the presence of both an 
rf trap potential and a strong magnetic field. An ion's 
momentum p is thus replaced by p— with the vector 

potential A satisfying B — \/ x A. The exact solution of 
the equations of motion contains quickly oscillating terms 
associated with the ion's micro-motion at the frequency 
fit of the rf trap drive. Averaging over the micro-motion 
will yield solutions characterized by an effective (pscudo- 
) potential [H,!!!]. 



i.e. a Hamiltonian for which x and y components are 
coupled with = eB/m. This coupling can easily be 
resolved by introducing new creation and annihilation 
operators 

a+ = —={ax + idy) und a_ = —={ax — iciy) , (16) 
V 2 V 2 

in terms of which the Hamiltonian reads 

("17) 

That is, analogously to the classical case, there are two 
decoupled modes with shifted frequency lOj. ± ^Wc, and 
the motion in z-direction is unaffected by the magnetic 
field. 
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V. DISCUSSION AND CONCLUSION 

In summary, we have proposed to encode quantum in- 
formation in nuclear spins of trapped atomic ions and 
considered the feasibihty of nuclear spin quantum infor- 
mation processing using trapped ^^^Yb+ ions in a lin- 
ear ion trap. Employing both nuclear and electron spins 
provides not only the combination of robust storage of 
quantum information with efficient quantum gating, but 
also a good way to suppress the undesired coupling be- 
tween electron spins. The discussion of possible methods 
to generate the required magnetic field indicates that this 
scheme is feasible with currently or near-future available 
ion-trap techniques. 

This scheme could also be applied to other candidate 
ions, such as '^^Ga+ d, O or ^Be 0. Since the hy- 
perfine coupling in those ions is much smaller than in 
171 Yb+, it is possible to satisfy the high-field approxi- 
mation using a lower magnetic field. On the other hand, 
because the nuclear spins of those ions are not 1/2 (i.e. 
43Ca+ with / =7/2 and ^Be with / =3/2), it would be 
more complicated to encode and manipulate qubits in 
nuclear spins. For example, as studied in the nu- 
clear spin flipping operation would be more complex and 
take a relatively longer time. Particularly, when x and 



y terms of the hyperfine coupling are considered in the 
Hamiltonian, a very high magnetic field is required to 
obtain the desired fidelity. Using the similar calculation 
for the CONTis gate for 7/2 nuclear spin ^^Ca+ in 0, 
we have found that 5T magnetic field are necessary to 
get an effective operation as good as that in the present 
paper for ^^^Yb"*". For ^Be+, at least 1 T magnetic field 
is needed to reach a good high-field approximation. 

In addition, without the magnetic field gradient, our 
scheme would still work using laser light for the elec- 
tronic spin operations using the Cirac-ZoUer model 
or M0lmer-S0rensen model [l^ . Due to involvement of 
nuclear spins, more qubits could be employed in the sys- 
tem with the same numbers of ions trapped compared to 
previous schemes using only electron spins. 
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